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We investigate classical and semiclassical aspects of 
codimension-two bifurcations of periodic orbits in Hamilto- 
nian systems. A classification of these bifurcations in au- 
tonomous systems with two degrees of freedom or time- 
periodic systems with one degree of freedom is presented. 
We derive uniform approximations to be used in semiclassical 
trace formulas and determine also certain global bifurcations 
in conjunction with Stokes transitions that become important 
in the ensuing difi'raction catastrophe integrals. 
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I. INTRODUCTION 

Periodic-orbit theory aims at the semiclassical eval- 
uation of energy levels of quantum systems and relates 
their spectral properties to periodic orbits of the corre- 
sponding classical system. For autonomous systems one 
considers the trace of the Green's function G{E) which 
determines also the density of states d{E). For periodi- 
cally driven systems the object of interest are the traces 
tr i^" of the stroboscopic time-evolution operator over n 
periods; they encode the so-called quasienergies of states 
that are stationary in the stroboscopic description. Both 
types of traces can be written as a sum of individual 
contributions of periodic orbits for chaotic (hyperbolic) 
systems Jnf-pj or a sum over rational tori for integrable 
motion pflqir 

Recent semiclassical studies based on periodic- 

orbit theory were devoted to the neighbourhood, in the 
space of control parameters, of classical bifurcations. 
These are instances in which periodic orbits coalesce and 
are the mechanism how orbits are born or disappear, or 
change their configuration when the energy or an exter- 
nal parameter is varied. Bifurcations are ubiquitous in 
systems with a mixed phase space and pave the path from 
integrable to chaotic motion. 

A collective treatment of the bifurcating orbits was 
found necessary, and even more the inclusion of predeces- 
sors of such orbits which live in complexified phase space 
and were termed ghosts. A collective contribution comes 
from an orbit cluster and can only far away from the bi- 
furcation be written as a sum of individual contributions. 
Both types of contributions are an additive term in the 
periodic-orbit expansion of the trace in question. 



The existing semiclassical (and most of the classical) 
studies focus on the generic bifurcations in the classifi- 
cation of Meyer and Bruno p^-p^ (see also |^). These 
are the bifurcations that are typically encountered when 
one has only a single parameter at hand to steer the 
system through parameter space, or, equivalently, when 
one investigates the periodic-orbit families in a given au- 
tonomous system as a function of energy. In general 
one assigns a codimension to each type of bifurcation 
by counting the number of parameters to be controlled 
in order to encounter it in a general setting. (The class 
of bifurcations of a given codimension is enlarged when 
symmetries are imposed on the system 
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generic bifurcations are accordingly the bifurcations of 
codimension one. In each of these bifurcations there is a 
central orbit of period n, surrounded by one or two satel- 
lite orbits of period nm. The cases m = 1, 2, 3, 4 are the 
tangent, period-doubling, period-tripling, and period- 
quadrupling bifurcations, respectively. There are two 
types of period-quadrupling bifurcations (island chain 
and touch-and-go), but only one for all other m. All 
period-m bifurcations with m > 5 follow the island chain 
pattern. 

In [0-^ transitional approximations for the collective 
contributions were derived that are only valid close to the 
(generic) bifurcation; far away from the bifurcation they 
give rise to individual contributions with the wrong am- 
plitudes. Refs. |10-12j give uniform approximations that 
are valid even far away from the bifurcation, where they 
asymptotically split into individual contributions with 
the correct amplitudes. 

The present work is devoted to bifurcations of codi- 
mension two in autonomous Hamiltonian systems with 
two degrees of freedom or time-periodic Hamiltonian sys- 
tems with one degree of freedom. The bifurcations are 
classified and their impact on semiclassical periodic-orbit 
theory is studied in detail. We derive uniform approxi- 
mations of collective contributions to the semiclassical 
traces and discuss certain global bifurcations in conjunc- 
tion with so-called Stokes transition. 

Bifurcations of codimension two manifest themselves 
in one-parameter studies in certain sequences of generic 
bifurcations. Sadovskii, Shaw, and Delos [^,^ found 
that such sequences can be explained by normal-form 
theory |p^-p6[|, but did not attempt a classification with 
respect to the codimension. The classical part of the 
present study is very much inspired by these works. 

It was indeed demonstrated in Ref. that bifurca- 
tions of codimension two are frequently felt semiclassi- 
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cally and necessitate a collective treatment even when 
one steers the system through control space with less 
than two parameters. This implies that collective con- 
tributions of this kind will constitute a basic ingredient 
in a semiclassical trace formula for systems with a mixed 
phase space. They indeed played an essential role in the 
semiclassical determination of the quasienergies of the 
kicked top |§. 

The paper is organized as follows. In section || we de- 
rive normal forms for the Hamiltonian that describe the 
bifurcations of codimension two. We find that they are 
organized by the multiplicity m in analogy to the situa- 
tion for codimension one. The corresponding sequences 
of codimension one bifurcations involve a tangent bifur- 
cation of period nm, followed by a period-m bifurcation 
that involves another orbit of period n. (The representa- 
tive case studied in corresponds to m = 3.) 

The normal forms and the sequences of codimension 
one bifurcations in the neighbourhood of the codimen- 
sion two point in control space (technically spoken, the 
unfolding of the normal forms) are discussed and illus- 
trated in sectio n ffl . 

With section |IV| we turn to aspects within semiclassi- 
cal periodic-orbit theory and present the starting point 
for the derivation of collective contributions of bifurcat- 
ing orbits, consisting of two-dimensional integrals over 
phase space or Poincare surfaces of section that involve 
the generating function S of the classical stroboscopic 
map. 

Uniform approximations of these contributions are de- 
rived in section They involve a phase function <i> 
and an amplitude function Normal forms of these 
are obtained from the corresponding Hamiltonians by 
non-canonical transformations and partial integrations. 
Investigating the influence from higher-order terms in 
the phase equips us with a sufficient number of coeffi- 
cients to guarantee the right stationary-phase limit of 
the expressions, which are therefore truly uniform. [The 
cases m = 1, 2 lead to standard diffraction integrals con- 
nected to the cusp and butterfly catastrophes, respec- 
tively. Among the large number of applications, the tran- 
sitional approximations have been investigated in con- 
nection to bifurcations of closed orbits in ||29| . Uniform 
approximations have not been derived there, however, 
and the canonical invariant determination of coefficients 
as well as Stokes transitions are also not discussed.! 



In section |V| we discuss certain global bifurcations 
that become important in the ensuing diffraction catas- 
trophe integrals. They give rise to Stokes transitions in 
which the contribution of a ghost satellite is switched 
on or off. The ghosts and transitions arise when the 
integrals are analyzed using the method of steepest de- 
scent. Stokes transitions have been investigated in the 
context of diffraction integrals and asymptotic expan- 
sions before. A uniform approximation for an isolated 
transition is given i n pof and has been applied for per- 
turbed cat maps in |31[, which is the only treatment of 
this phenomenon in semiclassics that we know of. The 



Stokes transitions investigated there, however, occur far 
away from any other bifurcation and can be regarded as 
isolated. A transition requires special treatment when it 
occurs in the immediate neighbourhood of a usual 'lo- 
cal' bifurcation, a situation that is often encountered in 
mixed systems. The uniform approximations and normal 
forms derived here can also be employed to describe the 
Stokes transition of a period-nrn ghost prior to a tangent 
bifurcation when the so-called 'dominant' orbit involved 
is real and of period n. The complete sequence of local 
and global bifurcations that we can handle consists of the 
period-m bifurcation at the central orbit and tangent bi- 
furcations of satellites, followed by Stokes transitions in 
which ghost satellites once more interact with the central 
orbit. 

We conclude and point out open questions in section 
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II. NORMAL FORMS OF THE HAMILTONIAN 
FOR BIFURCATIONS OF CODIMENSION TWO 

A. Objective 

The local bifurcations to be discussed are instances 
in which periodic orbits coalesce as parameters are var- 
ied. The types of bifurcations generically encountered 
in a given class of systems depends on the number of 
parameters varied, and the number of parameters typi- 
cally needed to be controlled in order to find a particular 
type is called its codimension. Here we investigate bi- 
furcations of codimension two in the class of periodically 
time-dependent Hamiltonian systems with one degree of 
freedom. In other words, we study families of Hamilto- 



H{q,p,t;e,a) 



(1) 



that depend on the two parameters e and a and obey 
H{t) = H{t+ 1), where the period is set to unity for con- 
venience. In general these systems have no time-reversal 
nor any geometric symmetry. The discussion directly car- 
ries over to autonomous systems with two degrees of free- 
dom since these can be reduced to one-parameter fami- 
lies of periodic systems with one degree of freedom by a 
standard procedure described e.g. in m. 



B. The bifurcation condition 

The periodic orbits show up as fixed points in iterations 
of the so-called stroboscopic map {q,p) {q',p') which 
is induced by the evolution over one period. This map is 
area preserving. Its linearized version 
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(corresponding to a 2 x 2 matrix) hence obeys det M = 1 . 
Orbits that appear for the first time in the rtoth iteration 
are said to be of primitive period hq . Such an orbit gives 
rise to uq fixed points in each n-step map with n — rrio, 
where r is an integer counting repetitions. 

The eigenvalues Ai,2 of the hnearized no-step map 
j\^("-o) are reciprocal to each other. A stable orbit has uni- 
modular eigenvalues and hence trM^"") = 2costi; with 
the real stability angle w. An orbit is instable if the 
eigenvalues are real. There are two cases depending on 
the sign of the eigenvalues, trM^""^ = ±2 cosh w', with 
the real and by convention positive instability exponent 
u;'. 

In general, an orbit bifurcates whenever the linearized 
n-step map M^"-* (with again n = rno) acts at the locus 
of the orbit in phase space in at least one direction as the 
identity map and hence obeys 



trikff™") =tr(Af("«))'' = 
or, equivalently, 

2 cos(27rZ/TO) 



trivr 



(3) 



(4) 



where the integers I, m are taken as relatively prime. 
This bifurcation condition implies a discrete m-fold ro- 
tational symmetry Cm in the flow pattern around the 
bifurcating orbit. For m > 2 the orbit in question is a 
'central' orbit on which 'satellite' orbits of primitive pe- 
riod nQin contract at the bifurcation. For m = 1 there are 
two possibilities, the orbit is either involved as a satellite 
in a bifurcation with an orbit of smaller primitive pe- 
riod, or it takes part in an isochronous bifurcation with 
other orbits of same period. Turning these observations 
around, there is always a central periodic orbit of smallest 
primitive period no among the bifurcating orbits which 
coalesces with satellites of period nom. For that reason 
m is called the multiplicity. 

The bifurcations of codimension one have been classi- 
fied by Meyer and Bruno [p^-p^ (see the introductory 
section). They constitute the building blocks of the sce- 
narios of higher codimension and will be illustrated to- 
gether with those of codimension two in the next section. 
Recall that for each m there is exactly one type with the 
exception of m = 4 which allows for two variants. 



C. Classification of normal forms 

The bifurcation condition (jj) is reflected by the Hamil- 
tonian flow around the central orbit; accordingly, the bi- 
furcations can be investigated by studying the general 



form of the Hamiltonian in the vicinity of this orbit. Fol- 
lowing Refs. [tP, ^ , ^ , ^ , we aim at the reduction of 
the general expressions to certain simple normal forms 
by suitable canonical transformations. In that way one 
can identify the parameters that govern the distance to 
the bifurcation. They can be chosen such that £ = 
brings us on a codimension one bifurcation and the codi- 
mension two scenario is encountered if in addition a = 0. 

For codimension one the construction (that is carried 
out for that case in detail in ^ and is recapitulated be- 
low) leads to the Birkhoff normal forms 



1 
2 
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> 5 



h^"'\q,p)-Ho 



eq + aq'^ + 

eq^ + aq^ + f 

el + a/3/2 3^ ^ 

el + aP{\ + cos 40) + hP{l - cos 40) 
aP + El'^.f biP + c/'"/^ cos m(j) 



zl 



(5) 



For m > 3 they have been expressed in canonical polar 
coordinates /, 6 with 



q = \f2I sin ( 



p — \/27cos( 



(6) 



The quantity is a constant. 

These expressions are autonomous and display the m- 
fold symmetry even globally. The periodic orbits are 
mapped to fixed points ^ = 0, ^ = and are thus de- 
termined as roots of polynomials in p and q. From Vieta's 
relations between these roots (or locations of satellites) 
and the coefficients of the polynomials it follows that or- 
bits collapse on the center as the lowest-order terms ~ e 
(and ~ a for codimension two) of the fixed-point equa- 
tions are steered to zero. 

To describe the codimension two variants we have 
to include higher-order terms to account for additional 
satellites that approach the center and obtain the ex- 
tended normal forms 
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a/3/2 cos 30 + 6/2 
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e/4 


a/2(l -1- cos 40) -1- 6/2(1 - cos 40) 






+cP{l + cos 40) 
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e/4 


a/2 + bl^/^ cos 50 
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e/4 


a/2 -1- 6/3 -1- c/3 cos 60 


> 7 


£/4 


a/2 + bil^ + c/™/2 COSTO0 



(7) 



The normal forms for to > 5 are the usual Birkhoff nor- 
mal forms; they will, however, be investigated not only 
for small e but also for small a. The expressions for m < 4 
go exactly one order beyond the Birkhoff normal forms. 



D. Derivation of normal forms 

In the derivation of the normal forms (||) and (^ the 
central orbit is placed into the origin of a local coordinate 
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system (g, p) by a time-dependent canonical transforma- 
tion. This is done in such a way that a Tayfor expansion 
in q and p yields 



Hiq,p,t)^Ho + -q^ + -p^ + Oi3) 



(8) 



with time-independent w and cr. The remainder in the 
expansion indicates third and higher orders in q and p. 
The new Hamiltonian will be ng-periodic but otherwise 
as general as the one we started with. Hence it suffices 
to study the case no — 1- 

One passes to a rotating coordinate system, where 
the angular frequency 27r//m is adopted to the motion 
around the center at the bifurcation, and examines the 
expansion of the Hamiltonian in a Taylor series in q and 
p as well as a Fourier series in t. Most terms in the 
expansion can be removed by canonical transformations 
(specified below) up to a certain order in p, q and the 
deviations of the parameters e, a that govern the dis- 
tance to the bifurcation. The expansion is carried out 
up to a certain degree, requiring that no bifurcating or- 
bits are added or qualitatively affected by the omitted 
higher-order terms. 



2. The cases m > 3 

For m > 3 the central orbit is stable close to the bi- 
furcation so that we can achieve u = cj in (|^), entailing 
trM = 2cosLi;. ft is then convenient to use canonical 
polar coordinates In the rotating coordinate sys- 

tem the leading-order term in I takes the form el with 
e = uj — 2tt— in agreement with the bifurcation condition 
. The expansion reads in detail 
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where m' runs from — fc,— fc + 2,...,A: since only such 
terms arise from expressions of type q^ pk-k ^ j^^^ 
assume that all t and dependence is already eliminated 
up to a certain order /'^/^. The majority of terms of this 
order are then removed by a canonical transformation 

to new coordinates J, 9 that is generated according to 

dG - 
di 



^ = J, ^ = 0, iJ' = iJ - ^ by the function 



1. The cases m = 1,2 

The linearized map AI has degenerate eigenvalues 1 
for m = 1 and —1 for m = 2. This entails that the cor- 
responding bifurcation scenarios of codimension one and 
two are essentially one-dimensional. The reason is that 
the linearized map is in these cases generically not di- 
agonalizablc: Otherwise one would have located the ma- 
trix M — ±1 in the three-dimensional manifold Sp{2, M) 
of real 2x2 matrices with detM — 1 Q. With only 
two parameters at one's disposal, however, one generi- 
cally finds only such matrices with tr M = ±2 for which 
the eigenspace is one-dimensional. The linearized map 
describes then a shear transformation and acts as the 
identity in only one direction. 

The reasoning can be put onto another footing by 
resolving the apparent paradox that we can cast the 
second-order terms into a diagonal form ojq'^ /2 + ap'^ /2 in 
the immediate neighbourhood of the bifurcation, which 
involves only two parameters. This diagonalization, how- 
ever, depends in a singular way on the parameters of the 
original expansion aiq^ +a2p^ +a3qp. (Note that all three 
terms are Ci- and C2-invariant.) Accordingly, three co- 
efficients flj are to be controlled in order to let the second 
order vanish. Such an argumentation will be used again 
for m — A. 

We choose the g-axis as the line on which orbits ap- 
proach each other at the bifurcation and invoke the split- 
ting lemma of catastrophe theory in order to simplify the 
expression in p-direction. With a suitable scaling trans- 
formation one can achieve \a\ = 1. The normal form for 
m = 2 is symmetric in q since in that case the coordinate 
system rotates by tt in each period. 



kl'm' 



(10) 
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2TTlm'/m - 2ttI' 



(11) 



After the transformation we switch back in our nota- 
tion from J, 9 to /, 4>. The coefficient Gwm' diverges 
at the bifurcation if the resonance condition I' = I— \s 

m 

met. This affects, for instance, all t- and ^-independent 
terms {V = m' = 0). The remaining 0-dependent terms 
are of type I^/"^ cos{nm(j) — 0fc„). Here n = 1,2,3,... 
is an integer since I and m are relatively prime, and 
k > nm as before. In the orders that appear in the 
normal forms the latter restriction admits only n — 1. 
The 0-dependent term of lowest order in / is generically 

of type I"^^^ cos (^m(f> + (pmij ■ A shift of </> eliminates 
the constant in the cosine. If the coefficient of this term 
is not small then one can get rid of constants (j)ki in 
higher orders A: > m by a transformation of the form 
<j) — 9 + Y^^'=i 9k'l'^ with suitably chosen coefficients. 



3. Further reduction for m = 4 

Additional considerations are needed for m — 4. The 
most general expression that goes one order beyond the 
Birkhoff normal form reads 
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=Ho + eI + aj2(l + cos4(/.) + 5/2(1 - cos 40) 
+ (c + d)I^ + (c - cos4(/) + e/^ sin 40 . (12) 

The codimension two bifurcation is approached for van- 
ishing e and a (or 6), while the other second-order co- 
efficient b (or a) is finite. Both cases are equivalent and 
mapped onto each other by a rotation about 7r/4. The 
normal form H^'^'^ has been written down for small a. 
Since b is finite we can eliminate two of the three third- 
order terms in (^2|) by a canonical transformation and 
achieve d = e = 0. [The corresponding generating func- 
tion is of the simple form G — 91 — {dP sin 49) / {8b) — 
eP/{8b) if corrections involving e and a are neglected; 
the complete form is slightly more complicated.] 

The fine-print in the derivation for to = 4 is that two 
orbits pretend to bifurcate also as we send a ^ b. We 
identify now these orbits and show that for codimension 
two they are actually ghosts (complex solutions of the 
fixed-point equations) at a finite distance to the center. 
For simplicity we set e = 0; this docs not affect the gen- 
eral line of reasoning. The satellites that concern us solve 
the fixed-point equation 



= -4/2[a-6+(c-d)/]sin40 = O (13) 



by / = /(°) = (5 - a)/(c - d). The other fixed-point 
equation yields 



cos-^ 



•-d 



(a - 6)2 



2^^CW. (14) 
a — 



(These satellites are related by the reflection symme- 
try ^ —0 and undergo a pitchfork bifurcation as 
|C*^°-'| = 1. The symmetry is broken if e 7^ 0.) For 
reasons similar to those put forward for to = 1,2, two pa- 
rameters have to be controlled in order achieve J'-^-' — 0, 
i.e., a = 6: For the given parameter combination there 
is no (j) dependence in the second order of /. But actu- 
ally there are two independent terms involving 0, namely 
cos 40 and sin 40 — one of them had been eliminated 
by a diagonalization that is again sensible only if the 
other has a non- vanishing coefficient. We already used 
two parameters, then, such that we must assume s and 
a « 6 to be finite. This gives | cos40| ~ (a — 6)"^ 3> 1 

and COS0 ~ (a - b)-^/^ = ©(/(o)"^^^), and the Carte- 
sian coordinates @ indeed remain finite: As announced 
this shows that the orbits that appeared to bifurcate are 
complex solutions (with real p and imaginary q) of the 
fixed-point equations and stay away from the center. 



III. LOCAL BIFURCATION SCENARIOS 

We discuss now in detail the bifurcations of codimen- 
sion two that are described by the normal forms given in 
the preceding section. In each case the location of the 



periodic points, given as the solutions of the fixed-point 
equations 



dH 
dq 







dH 

dp 







(15) 



are investigated as the parameters are varied. Sequences 
of codimension one bifurcations are encountered if only 
one parameter is varied close to a codimension two point 
1^,^. They are discussed here for fixed a and vari- 
able £ and are illustrated by contour plots of the normal 
forms. Unstable orbits appear there as saddles while sta- 
ble orbits correspond to maxima or minima. In all these 
sequences there is a period-m bifurcation at e = and a 
tangent bifurcation of satellites at the parameter combi- 
nations 
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> 7 



tangent bifurcation of satellites 



1 af_ 
4 62 

3 b , 

9 

32 b 

2 a_ 

128 

675 F" 
1 _af_ 

3 b±c 

3 b 



^ 00 L 



£ = — 
£ = — 



(16) 



Global bifurcations that are of particular interest in the 
context of uniform approximations are discussed in sec- 
tion fvl. 



A. Tangent bifurcations (m — 1) 

In a tangent bifurcation two orbits of the same primi- 
tive period coalesce. On one side of the bifurcation both 
orbits are ghosts, i.e., complex solutions of the fixed- 
point equations, and their coordinates and other char- 
acteristic quantities are related by complex conjugation. 
On the other side of the bifurcation both orbits are real, 
one of them being initially stable and the other unstable. 
The scenario is described by the normal form /i'^' which 
accounts for two periodic orbits ± at coordinates p± = 
and 



1 £ 

3 a 



(17) 



Often one encounters a third orbit of identical period in 
close neighbourhood (in phase space) to the bifurcating 
orbits. This orbit must be taken into account, for in- 
stance, to obtain a reasonable semiclassical approxima- 
tion. One has to work then with the extended normal 
form H'-^\ The fixed-point equation dH/dq = is a 
real cubic polynomial in q and has three solutions. The 
number of real solutions is determined by the sign of the 
discriminant 



D 



1 £ 

86 



1 £ /I a 
46 V4 6 



(18) 
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There are three real sohitions for D < and only one for 
D > which is then accompanied by two complex ones. 
Tangent bifurcations are encountered at Z? = 0, that is 
e^O or e = -a^/{Ab^). 

A sequence of these two tangent bifurcations is de- 
picted in Figure |l|. The codimension two bifurcation is 
obtained when s and a pass zero simultaneously. If this 
is done in such a way that the discriminant changes sign 
then the number of solutions changes from one to three 
in a pitchfork bifurcation. Such a bifurcation is even of 
codimension one if the system is time-reversal symmetric 
or has a reflection symmetry [|l6|-p0t . 

B. Period— doubling bifurcation (m = 2) and tangent 
bifurcation of satellites 

In a period-doubling bifurcation the central orbit 
changes its stability by absorbing or emitting a satel- 
lite of double period. In the Birkhoff normal form ft,*^^' 
the central orbit sits at coordinates Qq = po = and the 
satellite is represented by two fixed points 



.r = 0, .i=±f|^. (19) 

In the extended normal form ij'^' the central orbit lies 
again at qo ^ po = 0, but there are now two satellites ± 
with coordinates 

A tangent bifurcation of the satellites is encountered at 



but the condition ab < must be obeyed since other- 
wise both orbits are still ghosts with purely imaginary 
g-coordinates. A sequence of tangent bifurcation of the 
satellites and period-doubling bifurcation (with ab < 0) 
is shown in Figure |^. 

C. Period— tripling bifurcation (m = 3) and tangent 
bifurcation of satellites 

The situation for the period tripling is visualized in the 
sequence of contour plots in Figure 0. Initially, a stable 
periodic orbit of period one is surrounded by its stabil- 
ity island. At a certain value of the control parameter 
two satellites of triple period come into existence via a 
tangent bifurcation. Then the inner (unstable) satellite 
approaches the central orbit, collides with it in the pe- 
riod tripling, and finally re-emerges on the other side. 
This scenario has been investigated, for instance, in the 
diamagnetic Kepler problem and for the kicked top 

i3- 



The Birkhoff normal form /i^^^ describes the central 
orbit at / = and the unstable satellite that is in- 
volved in the tripling. The (/i-coordinate of the satellite 
obeys dh^^^ /dcj) = — 3a/^/^sin30 = 0. Since a three- 
fold symmetry is implied by this equation it suffices to 
consider the second equation dh'-^^ /dl = on the p- 
axis after switching back to the coordinates p, q, yielding 
ep + -^CLP^ = 0. Hence the canonical radial coordinate 
of the satellite is / p^ 12 = Ae^K^a^). 

In the extended normal form H^^^ the ^-coordinates 
of the satellites once more obey -3a/3/2sin30 = 0. On 
the p-axis they satisfy now ep + "^"p^ -f- bp^ = 0. This 
equation has three solutions, 

3 a r%a? £ 

One in fact sees that the inclusion of the next-order 
term implies the existence of a further satellite. At 
e — 9a^/ (32&) the satellites undergo a tangent bifurcation 
and for e/b > 9a^/(326^) both satellites are ghosts. For 
Q<e/b< 9aV(3262) both satellites are on the same side 
of the central orbit, while after the period tripling (e = 0) 
they lie opposite to each other. In the limit e/b — oo 
the satellites form a broken torus, well separated from 
the central orbit. When a second parameter is varied to 
achieve a = e = 0, both satellites are contracted onto the 
central orbit in the codimension two bifurcation. 



D. Period— quadrupling bifurcation (m = 4) and 
tangent bifurcations of satellites 

There are two variants of the period-quadrupling bi- 
furcation depending on the magnitude of the coefficients 
a and b in the normal form h^'^\ In both cases there 
are two satellites of quadruple period involved that lie 
at sin 40 = and are distinguished by the quantity 
cos 40 — ±1 = a. Their radial distance is given by 
/(^=i) ^ -£/(4a) and /("^^-i) = -e/(4fe). In the touch- 
and-go case sign a = — signfo an unstable satellite be- 
comes a ghost while in turn a ghost solution becomes real 
and emerges from the central orbit. In the island-chain 
scenario sign a = sign 6 there are two ghost satellites on 
one side of the bifurcation and two real satellites on the 
other, one of them being stable and the other unstable. 

The next-order terms in the extended normal form 
ijf") [eq. (|l|)] involve three new parameters c, d and 
e and give rise to six satellites. For e = there are 
two satellites on the lines cos 40 = 1 and two on the 
lines cos 40 = — 1 as well as the two satellites discussed 
in the derivation of the normal form H^^\ There is a 
tangent bifurcations at 2a^ = 3ec where the satellites on 
cos 40 = 1 coalesce, and another one at 2b^ — 3ed that 
involves the satellites on cos 40 = — 1. A great variety 
of possible configurations of all six satellites exists. Here 
we are, however, only concerned with the codimension 
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two bifurcation, described by ff^"*) and encountered for 
£ = a = 0. It involves only three satellites, those on the 



lines COS' 



1 with radial coordinates 



(1) 



la ^ la? 1 e 



(22) 



and that satellite with cos 4^ = — 1 which is closer to 
the center and lies with 7?^^^ at = — £/(4&). Com- 

pared to the situation described by H^'^^ the second satel- 
lite on the line cos 40 = — 1 is shifted to infinity; the 
two satellites at / = (6 — a)/d have now angular coordi- 
nates cos 40 w — 5 and are therefore ghosts. (Certainly 
they eventually may become real at finite values of e and 
a, far away from the codimension two bifurcation and 
therefore out of the scope of the present work.) A tan- 
gent bifurcation is met at e = 2a^/(3c) provided that 
/ — — a/(3c) > 0, since the Cartesian coordinates (||) 
are otherwise imaginary. Sequences of a tangent bifur- 
cation at positive / and the two variants of quadrupling 
bifurcations are shown in Figure |4[ 



E. Period— m bifurcation with m > 5 and tangent 
bifurcations of satellites 

The codimension one bifurcations for m > 5 follow 
the island-chain pattern already encountered for m = 4: 
There are two satellites that are ghosts on one side of 
the bifurcation and real on the other, one of them being 
stable and the other unstable. The stable and unstable 
periodic points form a chain similar to the broken rational 
tori that appear in almost integrable systems. Indeed, 
the 0-dependent terms in the normal forms are of the 
form of a small perturbation in that situation. 

For m > 5 the usual Birkhoff normal forms describe 
even the codimension two bifurcations: in addition to 
the orbits participating in the period-m bifurcation of 
codimension one they also account for the satellites that 
are involved in the subsequent tangent bifurcations. In 
the case m = 5 one obtains three satellites at sin 50 = 
that satisfy on the p-axis 



ap 



4%/2 



6^=^ = 



As for TO = 1 it is the discriminant 

I A y \ ^ 9. p / R n \ ^ 

D 



2e 
56 



8 a 
15 6 



(23) 



(24) 



of the equation that governs the number of real solutions. 
There is the period-5 bifurcation at e = and a tangent 
bifurcation at e = — 128a'^/(6756^). That sequence is 
depicted in Figure ^. 

In the case m — 6 one finds four satellites, two on each 
of the lines cos 60 = ±1 = cr at 



(25) 




Tangent bifurcations take place at independent parame- 
ter combinations e = a'^ /[3{b + ac)] (provided that the I- 
coordinate is not negative) . A sequence with two tangent 
bifurcations at positive values of / is shown in Figure |^. 
All four satellites approach the center in the codimension 
two bifurcation as e and a are sent to zero. 

For TO > 7 there are even more satellites in the Birkhoff 
normal form than the four that participate in the codi- 
mension two scenario. In first order, the relevant satel- 
lites lie on the lines cos to0 = ±1 = cr at a radial distance 



la2 
9 62 



l£ 

3 6 ' 



(26) 



which is independent of a. The 0-dependent term in- 
duces a small correction of order (t/'"/2-2 Before the 
tangent bifurcations, which are encountered at almost 
identical values e « 0^/(36), both satellite pairs have 
complex /. After the bifurcation both the inner as well 
as the outer orbits form island chains that are visible in 
phase space if the /-coordinate is positive. The corre- 
sponding sequence is shown in Figure ^. 



IV. PERIODIC ORBITS IN SEMICLASSICAL 
APPROXIMATIONS 

Individual contributions of periodic orbits to semiclas- 
sical expressions of traces lose their validity close to bifur- 
cations. A collective treatment of the involved orbits is 
then necessary. Here we discuss general aspects [0-^|j8| in 
which uniform approximations that provide the required 
regularization are integrated in the next section. 

Gutzwiller's trace formula relates the trace of the re- 
tarded Green's function 



G{E) 



I 



E + iO+ - H 



(27) 



of an autonomous system with Hamiltonian H to the 
properties of the periodic orbits in the classical system 
and thus provides a semiclassical expression of the den- 
sity of states 



d{E) = --Im tTG{E) 

TT 



(28) 



A suitable starting point for the derivation of periodic- 
orbit contributions is given by an integral over Poincare 
surfaces of section 17 (see e. g. |lo[]). 



n=l 

X exp 



dqndpo 



1/2 



^[S^"\qn,Po;E) - QnPo] ~ i^l' 



(29) 



This involves the generating function S^'^\qmPQ]E) of 
the nth iterate of the Poincare map (goiPo) ~^ {q-mPn) 
with 
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90 , 



BE 



T, (30) 



where T denotes the time elapsed along the trajectory. 
The integration domain Q might be disjunct and is also 
used to account for the different sheets of the generat- 
ing function, i.e., the multivaluedness of S^'^\qn,Po', E). 
Finally, v is the Morse index. 

A similar expression can also be obtained for periodi- 
cally driven systems that are described stroboscopically 
by a time-evolution operator F. For convenience we set 
the stroboscopic period to unity. F represents in general 
the unitary operator of a quantum map with a classical 
limit. The eigenstates of F are stroboscopically station- 
ary, and the phases of the unimodular eigenvalues are 
called quasienergies. The quasienergy spectrum is en- 
coded in the traces trF". With the Van Vleck propaga- 
tor one obtains the expression 



dg„ dpo 



2T:h 



Q2g(r 



Indpo 



1/2 



X exp 



^[S'("^(q„,Po) - qnPo] - 



(31) 



where S'^"^(g„,po) is now the generating function of the 
n-step map with 



dqn 



= Pn 



dpo 



= qo 



(32) 



and the integration domain lies in phase space [again, 
accounting also for the multivaluedness of 5'("^((7„,po)]- 
We consider contributions to trF" or ih tr G{E), the 
factor being introduced to facilitate a parallel investiga- 
tion of both cases. Let us examine the contribution of an 
arbitrarily chosen region O', 



2'Kh 



= TT^ / dg' dp *(g',p) exp 



(33) 



with the notations q' = g„, p = po- For the contributions 
to ihtiG{E) the phase function $ and the amplitude 
function ^ are 



'^{q',p)^S^-\q',p;E)-q'p 

1/2 

^{q',p) 



1 as-^") 

n dE 



dq'dp 

For contributions to tr F" they read 

<P{q',p)^S^-Hq',p)~q'p 
1/2 

*(g',p) = 



dq'dp 



(34) 
(35) 

(36) 
(37) 



The main contributions to the integral ( |33| ) arise from 
regions Q' around the stationary points. These are given 



by ^ = and ^ = which corresponds to 



dq' 



dp 



dpo 



= Po ■ 



(38) 



The solutions are the periodic points of a given energy 
in the autonomous case or given stroboscopic period in 
the driven case. These regions come naturally into focus 
when one uses the steepest-descent method to find the 
leading-order term of an asymptotic expansion of the in- 
tegral in H. To achieve this goal the integration variables 
are complexified and the initial contour is deformed such 
that the maxima lie at the solutions of (|3^). From the 
maxima one follows paths of steepest descent of the inte- 
grand. The new contour has to originate from the original 
one by a continuous deformation without crossing singu- 
larities. In order to construct a contour that connects the 
original integration boundaries one automatically visits 
also some complex 'ghost' solutions of (^8|). Only the 
ghosts that lie on the deformed contour are relevant. We 
come back to this point in the discussion of global bifur- 
cations in section 



The size of the region fl' of almost stationary phase 
depends on h and is determined by the condition that the 
typical variation of the phase ^{q\p)/h over that region 
is of order 27r. The region can contain several stationary 
points, i.e., phase space points of a number of periodic 
orbits, and each orbit may contribute with several points. 

The region shrinks if h is sent to zero while all param- 
eters are fixed at general values (codimension zero), and 
finally splits into regions containing only a single periodic 
point. Each of these gives rise to a usual stationary-phase 
(sp) contribution 



(sp) 



To exp [j-SM _ q^-^ 



no 



tr M("o'')p/2 



(39) 



There are no such contributions from each orbit. More- 
over, repetitions of an orbit are regarded as independent 
here. Four canonical invariant characteristic quantities 
of periodic points enter this semiclassical expression, the 
primitive period Tq, the action 5*^"^ the stability factor 
trM("\ and the Maslov index fi. They also determine 
the uniform approximations to be derived in the next 
section. 

In the periodically driven case. To = no is the primitive 
stroboscopic period that was introduced in section ||. 
(Recall that we set the stroboscopic period to unity.) The 
action S is given by the value of Hamilton's principal 
function, S — J{pdq — Hdt), or, equivalently, by the 
value of the phase 



S^^\q,p) 



qp 



(40) 



at the periodic point. The linearized n-stcp map M*-"-* 
was likewise introduced in section |^ It is connected to 
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the second derivatives of S**^"^ and involved in the expres- 
sion through its trace 



trAf(") = 



1 + 



\ dq' dp J 



dq''^ 



dq' dp 



(41) 



FinaUy, there is the Maslov index fjL — v— ^sign G where 
V is the Morse index and signG denotes the difference 
in the number of positive and negative eigenvalues of the 
matrix 



G = 



9^$ d'^-i' 

dq'^ dq' dp 

d^<s> a^* 

dq' dp dp^ 



(42) 



involving second derivatives of $ = S''^"^ — g'p. 

The four quantities are almost of the same form in the 
autonomous case. The primitive period Tq is the small- 
est period T [see eq. (|o|)] after which one comes back to 
the initial point when starting somewhere on the trajec- 
tory of the periodic orbit. The orbit shows up as a fixed 
point in all n-step maps on the surface of section where, 
as before, n — rno, and r = T/Tq is an integer count- 
ing repetitions. The quantity S is again given by ( |40| ) 
but corresponds now to the reduced action S = §pdq. 
The linearized map on the surface of section and the sta- 
bility factor are connected to the generating function in 
the same way as before, and the Maslov index is again 
connected to the Morse index and the matrix (^) by 



isignG. 



V. UNIFORM APPROXIMATIONS 

A. Breakdown of stationary phase near bifurcations 

From the bifurcation condition (^) it follows that the 
individual contribution (^9|) of an orbit blows up close to 
a bifurcation and even diverges right at trM = 2. The 
assumption under which a stationary-phase approxima- 
tion is reasonable is then no longer fulfilled: It does not 
suffice to expand the generating function up to second 
order around the trajectory. 



$(<Z',P) -5*0 



UJ 



-p2 + 0(3) 



(43) 



since the coalescing orbits cannot be separated even in 
the limit h 0. In the two-parameter family of Hamilto- 
nians (^) the stationary-phase approximation is in dan- 
ger close to bifurcations of codimension one and two. In 
practical applications with small but finite values of h, 
one typically even visits regions in which bifurcations of 
higher codimension are felt. Along the same line, bi- 
furcations of codimension two are even felt when only a 
single parameter (frequently, none at all) is varied; this 
observation is indeed our principal incentive. 



In this section we derive normal forms for the phase 
function $ and the amplitude function ^E" in (|3^) that 
supersede the quadratic form ( ^ ) and yield regular ex- 
pressions as a substitute for the stationary-phase result 
( ^ ) in regions affected by bifurcations of codimension 
two. 



B. Normal forms for phase and amplitude 

To each normal form of the Hamiltonian there is a cor- 
responding expression for the generating function which 
carries over to the phase function This function has 
as many stationary points as the Hamiltonian. The sim- 
plest functional form that can be achieved is identical to 
the normal form of the Hamiltonian, but with altered co- 
efficients. (Observe, however, that although this form of 
^(q',^), when expressed by canonical polar coordinates, 
obeys again rotational symmetries Cm, this is no longer 
the case for the map generated by it.) To illustrate this 
identification we note that the normal forms effectively 
describe the integrable dynamics of an autonomous sys- 
tem with one degree of freedom. In action-angle variables 
J, tp the evolution over a time interval of duration one is 
generated by 



5( J, V') = 5*0 + + 27rn) J - H{ J) , 



(44) 



such that J = J' and -0 = -0' — w( J) mod 2tt with the 
torus frequency a;(J) = dH/dJ. 

The action variable J is quantized since the angle ip is 
only given modulo 27r; this gives rise to branches of the 
generating function that are here enumerated by n. We 
circumvent this obstacle by considering the map in Carte- 
sian coordinates, generated by S{q' ,p), and appealing to 
the canonical invariance of the leading-order term in the 
?i-expansion that we are looking for. [This does not affect 
the option to shift finally to canonical polar coordinates 
in the integral expression (|3^).] The transformation to 
these variables yields new coefficients: for instance, while 
the second derivatives of H around a stable orbit involve 
the bare stability angle w, we are led for the generat- 
ing function to the relation ( ^l|) with trM = 2cosa-'. 
The coefficients in <& are therefore related, but not iden- 
tical to those in the normal forms of the Hamiltonian, 
although we will not refiect this in a change of notation; 
however, minus signs will be introduced following a con- 
vention that is motivated by eq. (^). 

In section || we observed that the normal forms de- 
scribe in some cases additional orbits that do not partic- 
ipate in the bifurcations of codimension two. In addition 
we will see that there are in general not enough coeffi- 
cients to yield independent actions and semiclassical am- 
plitudes in the stationary-phase approximation. To over- 
come these restrictions one has to consider the influence 
of higher-order terms in the normal forms on the classical 
properties of the orbits. These terms can be eliminated 
in the region of almost stationary phase by non-canonical 
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transformations. Such transformations can also be used 
to get rid of the additional non-bifurcating orbits de- 
scribed by the original normal form. The Jacobian of the 
transformation enters the amplitude function which 
can be simplified further by par tial integrations. The 
procedure is sketched in section VP . It results in the 
following normal forms for $ and 



m 


p) - 


So 




1 


—eq''' — aq''^ 


- hq'" - f 




2 


—eq — aq* 


- bq'^ - f p2 




3 


-el - a/3/2 


cos3(/>-6/2 




4 


-el -aP{l 


+ cos 4(/)) - hi 


^{1 - cos 40) 




-c/3(l 


f cos 40) 




5 


-el - aP - 


bPP cos 50 ~ 


ecP^^ cos 50 


6 


-el - aP - 


bp — cP cos 6 


)0 — edP cos 60 


> 7, odd 


-el - aP - 


bP-p/^{cI 


+ ed) cos TO0 


> 8, even 


-el - aP - 


bP-P{cI + 


ed) cos TO0 



(45) 



m 


<5{™)(q'^p) 


1 


1 + aq' + Pq''^ 


2 


1 + aq'^ + 


3 


l + al + (3P 


4 


1 + aI + jSP + 7/2 cos 40 + SP 


5 


1 + (a + /3/)/i/2 cos 50 + 7/ + 5P 


6 


1 + (a + 131)1 cos 60 + 7/ + 5P + e,P 


> 7, odd 


1 + (a + /3/)/i/2 cos m0 + 7/ + 6P + S,P 


> 8, even 


1 + (a + 131)1 cos m(j) + ^I + SP + ^P 



(46) 

In we demand \a\ = 1. The normal forms for m > 3 

are expressed in canonical polar coordinates /, defined 
in equation (|^). Note that terms show up in the expres- 
sions for m > 5 that cannot be expressed 'perturbatively' 
as in Cartesian coordinates (^. The substitutions 
m0 = 2'0 for m even and m0 = tp for m odd provide a 
potentially useful regularization. 

Collective contributions of bifurcating periodic orbits 
are obtained by inserting the normal forms for $ and ^' 
into (H). 

Numerically useful expressions in terms of Taylor series 
can be found for some of the integrals in |2^] {m — 1, 2) 
and 1^ (m = 3). The integrals for m = 1, 2 can also be 
easily evaluated by the method of steepest descent (cf. 
section VI) since they are essentially one-dimensional. 



For m > 3, however, a two-dimensional steepest-descent 
manifold might be quite difficult to construct. It is per- 
haps more convenient to deform only the / coordinate 
into the complex, yielding a simple steepest-descent con- 
tour for each fixed, real 0, and then to perform the 0- 
integral of finite range. 



C. Determination of coefficients 

The coefficients in the normal forms have to be ex- 
pressed by the classical quantities S, trM of the orbits 
in order to obtain a contribution that is invariant under 
canonical transformations 



la 1 



35| 



This is achieved 

by examination of the stationary-phase limit /i — *■ of 
(33) while all other coefficients are fixed. With definition 
(42), each stationary point gives rise to a contribution 



exp 



— 5> — I— { ly signG 

h 2 \ 2 ^ 



, (47) 



which determines the parameters by comparison to (|3£ 
The collective contributions are then not only applica- 
ble in the immediate neighbourhood of the bifurcation, 
but also far away, where they split into a sum of isolated 
contributions of Gutzwiller type ( p9| ) with correct am- 
plitudes and phases. Consequently they constitute uni- 
form approximations. (Transitional approximations of 
the type mentioned in the introduction are obtained if 
one uses ^ = 1 instead.) 

In detail, the properties of the central orbit determine 
e and So, since the stationary point in the origin gives 
C^-P) = exp [iSo/h \{(T + signe))] / .J\2l\ for 

TO = 1,2 (recah that |cr| = 1) and C^'p^ = 
exp \j.So/h — i-^iv + signe)] /\e\ for m > 3. The remain- 
ing coefficients of the phase function are uniquely deter- 
mined by the actions of the satellites. It turns out to 
be helpful to use an ansatz where the coefficients are ex- 
pressed by scaled positions on a radial line connecting 
the satellites with the central orbit. For even to with two 
real satellites on such a line, for instance, we put them 
on scaled positions dX xi — ±1, X2 = ±?/, corresponding 
to 



d$ 

dx 



(48) 



integrate to obtain <^{x) and determine y from the (scale 
invariant) ratio 



Si — So 



l-3y2 



^2 - ^0 y\y^ - 3) 



(49) 



Without restriction we can demand < y < 1; then 
there is exactly one solution. The factor A follows from 
the absolute value of 5*1 — S'o, and the scale of x is fixed 
by knowledge of e. In the case of complex satellites they 
are placed at a; = ±1 ± iy, i. e., $ = Ax'^[x^ + 3x^{y^ — 
1) -I- 3(?/^ -I- 1)^], and y is obtained from 



RcSi-So 1 + 9y2 _ - 



Im 5*1 



16?;3 



(50) 



There is a solution with \y\ < 1 and another one with 
\y\ > 1. The right choice takes into consideration whether 
the ghost with ImS" > lies on the steepest-descent 
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contour connecting the integration boundaries or not; see 
section 



V] 



The approach presented here to obtain the coefficients 
of $ works also for the other normal forms. Moreover, 
the stationary-phase result is a linear combination of the 
coefficients of the amplitude function which are there- 
fore easily obtained by comparison with the semiclassical 
amplitudes in (^). For m > 5 there is a symmetry- 
related pair of 'spurious' ghost satellites (analogous to 
those already discussed for m = 4), which are negligi- 
ble not only since they do not lie on a steepest-descent 
contour (see again section Vl), but also because a 'magic' 
cancellation in course of the derivation (see below) entails 

= at their positions. Fortunately, exactly one extra 
term in ^' shows up in these cases which can be used 
to achieve this suppression. It seems reasonable that one 
uses this approach also in the case m = 4, where an extra 
coefficient is also at one's disposal. 



D. Derivation of normal forms 

We perform now the reduction of the original normal 
forms for $ (identical in appearance to those of H) to 
the forms listed above and obtain in parallel the expres- 
sions of 4*. The remainder consists of higher orders in 
/ as well as a, e. Once more we invoke Vieta's relations 
and regard the coefficients as certain orders of the typical 
distance of the satellites to the central orbit. The coeffi- 
cient £ is usually the product of such distances, while a 
is a sum. We use that the original normal forms entail 
^' = const plus corrections of the form q''^p''. In most 
cases the precise form of coefficients in the transforma- 
tions is of no particular interest and therefore not given; 
the expressions are easily obtained, for instance, with the 
assistance of symbolic mathematical programs. Only the 
feasibility of reduction counts. Finally all coefficients in 
^ and $ are not to be determined from explicit expan- 
sions of Hamiltonians or generating functions but rather 
from the actions and stability properties of the orbits as 
explained above; otherwise no canonically invariant re- 
sult would be obtained. 

For m = 1, 2, 3 all the periodic points described by the 
original normal forms are involved in the bifurcation. For 
m ~ 1 there is at least one real orbit which we place into 
the origin by a shift of q'. This results in the normal 
form <i>(^). The three coefficients e, a, and b as well as 
the value So = 4'^^'' (0, 0) are fixed by the three actions 
of the orbits and one of the stability factors. From the 
normal form follows ^E* = 1 , the other stability factors are 
therefore not yet independent. At a little distance to the 
bifurcation, however, higher-order terms in the Hamilto- 
nian or the phase function act as a perturbation, and the 
implied relations between the classical quantities of the 
orbits are no longer valid. For m = 1 we consider terms 
of type cq'^ -\- dq'^ . Terms involving p effectively do not 
alter the final expression and can be discarded by ap- 



pealing to the splitting lemma of catastrophe theory. In 
the region of almost stationary phase the higher-order 
terms act as a perturbation and can be eliminated by 
substituting q' = Q + AQ^ + BQ^ with suitably chosen 
coefficients. The Jacobian of the transformation involves 
Aq' = AQil + 2AQ + 2>BQ^) and gives the normal form 
^'"^\q' tP) announced above. The two additional coeffi- 
cients are determined by the remaining stability factors 
tr M of the satellites. Corrections to $ of even higher 
order would carry over to higher-order terms in Vf. They 
involve additional coefficients and on first sight allow for 
ambiguities, but can be eliminated by successive partial 
integrations. The term of highest order is written as 



;-3 



dg 



terms of order I — \, I 



(51) 



The partial integration of the first term gives an order 
hQ^~'^ and a boundary contribution that vanishes for I > 
4. In the course of this procedure the constant 1 and 
the coefficients a and (3 in ^^^^ acquire next-to-leading 
order corrections in h that can be discarded and reflect 
canonically non-invariant properties of the orbits. 

For TO = 2 the higher-order terms have to obey the 
reflection symmetry q' —q' and are on the g'-axis 
of the form dq'^ + eq'^^ . They are eliminated by q^ = 

+ AQ'^ + BQ^. The Jacobian involves only even orders 
of g. Terms of order I > 6 can be eliminated in ^! by 
writing 



g' ~g 



/-5 



dg 



-I- terms of order / - 2, Z - 4 (52) 



and performing a partial integration of the first term that 
yields a term ~ hQ^^^ . This reduces the amplitude func- 
tion to its normal form 

In the case to = 3 we have to consider the influence of 
the terms c/^/^ cos?,(l)-\-dP . We can safely use e = ©(/*) 
and a — 0{I*^^^) as upper bounds in orders of the typ- 
ical distance /* of the satellites to the center, which is 
in turn connected with the size of the region of almost 
stationary phase. Note that this does not impose a re- 
striction on the relative order of these parameters as long 
as they are small enough. Having this in mind, it is 
easy to see that the elimination of the extra term can be 
achieved by a transformation I = J + AJ^^^ cos 3^^ and 
(f> = ip + BJ^/^ sin 3-0. The corresponding coefficients A 
and B give rise to a cancellation of the order J^/^ in the 
Jacobian. This order also does not show up in the orig- 
inal since it can expressed by terms of the form q'^p^ 
and obeys a three-fold symmetry. Hence, the order /^/^ 
is absent in Higher orders are again eliminated by 

partial integrations, carried out now with respect to /. 
The integration over suppresses terms that are odd in 
(/>, such as sin 30. Note that one could work alternatively 
with the form ^(■^^ = 1 + al + cos 30; the equiva- 

lence to '^'^^^ is again worked out by a partial integration. 

Only three of the six satellites described by H'^^^ ( p^ ) 
are involved in the codimension two bifurcation with 
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m = 4. The normal form <I>(^) given above has been 
given a deeper foundation in the discussion of the corre- 
sponding Hamiltonian H^^\ The outer satelhte on the 
hues cos4V' = — 1 is shifted to infinity and is infinitely 
unstable, |trM| — oo. The two satellites at radial dis- 
tance / = {b — a)/c have cos4V' ~ —5. They are con- 
sequently ghosts with real actions and do not contribute 
in the stationary-phase limit (for a deeper foundation 
see section Vl); moreover, they are quite far away from 
the bifurcating orbits. We assume for that reason that 
their influence is negligible. This leaves us with a nor- 
mal form that is again completely determined by the ac- 
tions and one stability factor. We assume that the three 
non-bifurcating satellites remain negligible even under 
the influence of higher-order terms. The corresponding 
amplitude normal form ^'^'^^ is irreducible under further 
partial integrations, but has one more coefficient than 
needed to account for independent stability factors of the 
satellites. In analogy to the situation to be discussed for 
m > 5 we can use this coefficient to yield ^ = at the 
position of the unwanted ghosts in favor of an additional 
suppression. 

For m = 5 there are three bifurcating satellites but not 
enough coefficients in the original normal form to account 
for independent actions. Observe that a scaling transfor- 
mation / = AJ does not affect the values of the phase 
function at the stationary points. Only the two combi- 
nations e^/a and enter these values which have to 
match the three actions of the satellites. (The coefficients 
are determined uniquely if one takes the stability factors 
into consideration.) Independent actions are admissible 
after allowing for higher-order terms which are once more 
removed by a transformation in order to yield no spuri- 
ous additional stationary points. The next-order term, 
conveniently expressed as -|-56c/2/^, can be eliminated 
by a transformation 



/ = J + cos Si/' 



(53) 
(54) 



which is similar to the one for to = 3, but now the order 
J^/^ survives in the Jacobian, 



* = 1 - cJi/^ cos 5?A-K ©(J) 



(55) 



The transformation gives rise to the term —ecfil"^ cos hi\) 
in $ and provides us with the additional coefficient c. 
[For illustrative purposes we wish to mentioned here that 
the coefficient ec in $ can be treated as C'(J*'^^^), where 
J* gives the order of the distance of the satellites to 
the central orbit. This order is related to the coeffi- 
cients e and a by application of Vieta's relations to the 
stationary-point equation d^/dl = 0.] The next-order 
corrections to $ give even one coefficient more in than 
necessary for independent stability factors of the bifur- 
cating orbits. There is, however, an extra pair of satellites 
at / = /(o) = -ec/h that even approaches the center as 
£ ^ 0. The angular coordinate, however, obtains a large 



imaginary part since it obeys (/^°^)^/^ccos50 = 1, where 
a term of order a has been dropped. Indeed this yields 
in leading order ( |55| ) a vanishing \1/ and encourages us to 
use the extra coefficient to accomplish suppression of the 
unwanted ghosts. 

A similar situation is encountered for m — Q: Only 
three independent actions of four bifurcating satellites 
can be modeled with the original normal form, but 
higher-order terms give rise to corrections that lead to 
the given normal forms of phase and amplitude. An ex- 
tra coefficient is again present to suppress the ghost pair 
at / = — ed/c, and the Jacobian of the transformation 
turns out to be once more in favor of such a strategy. 

Enough coefficients for the actions are present in the 
original normal form of $ for to > 7, but the expres- 
sion yields more stationary points than desired. The 4>~ 
dependent terms of highest order in / can be eliminated 
in favor of terms of lower order by successive substitu- 
tions I ^ J + AJ^ cos nKj). This procedure can be carried 
out in parallel to substitutions / = J + BJ^ that aim at 
the elimination of (/(-independent terms. The coefficients 
of the remaining terms reflect the values at the station- 
ary point up to the order that yields them as independent 
from each other and allows also for independent stability 
factors through the expression for In the derivation 
one encounters again a Jacobian that vanishes at the lo- 
cation of the unwanted ghost pair at / = —ed/c. We 
should note that the coefficient c is here of order a, or, 
equivalently, e^/^. 



VI. STOKES TRANSITIONS 



A. Preliminary remarks 



We mentioned in section IV that a steepest-descent 
contour has also maxima in complexified phase space 
which correspond to ghost solutions of the stationary- 
point equations. A helpful rule in that respect is that 
only ghosts with Im $ > can lie on the steepest-descent 
contour. Moreover, satellites that disappear in period-TO 
bifurcations of codimension one with even m are after- 
wards 'self-conjugated' ghosts, that is, map onto them- 
selves under complex conjugation, and have therefore real 
classical quantities (in canonical polar coordinates they 
have real (j) and / < 0); accordingly, they do not con- 
tribute. The ghosts immediately beyond a period-TO bi- 
furcation of codimension one with odd m > 5 have almost 
real action, a small imaginary part only being introduced 
from higher orders, and also do not contribute. For the 
remaining ghosts, however, it cannot be avoided to con- 
struct the contour in order to find out whether they are 
relevant or not (although the majority of relevant ghosts 
will be close to reality, i.e., about to bifurcate). 

A steepest-descent contour consists of different sheets. 
The phase Re $ in the exponent of the integrand is con- 
stant on each sheet and thus given by its value at the 



12 



periodic point, that is, the real part of its action. For 
general combinations of the control parameters (codi- 
mension zero) there will be only one orbit lying on each 
sheet, though it is possible that it does so with more 
than one of the points along its trajectory. Imagine now 
that for one combination of the parameters a ghost lies 
on the steepest-descent contour while for another one 
it does not. The ghost is denoted by + in the follow- 
ing and sometimes called 'subdominant' orbit. Its com- 
plex conjugate partner is denoted by — . Somewhere on 
a path connecting both parameter combinations the con- 
tour changes its form qualitatively in a so-called Stokes 
transition: The sheet of the ghost merges with the sheet 
of another orbit, which is called the 'dominant orbit' and 
denoted by 0. On both sides of the transition the sheets 
of the two orbits will connect different zeroes of the in- 
tegrand. In the transition the contour changes in such a 
way that the sheet of the subdominant orbit is no longer 
needed to connect the original integration boundaries. 

A necessary condition that the sheets of two orbits 
merge is that the real part of their actions are identical. 



Re So = Re 



(56) 



In general, this condition is not sufficient since both 
sheets could be separated by others. (Investigating ener- 
gies of real orbits, the analogous condition Hq = H+ 
is necessary to find a heteroclinic orbit; again, addi- 
tional insight is needed to decide whether an equipoten- 
tial contour joining both orbits indeed exists. These two 
global bifurcation types are therefore intimately related 
and constitute, in the language of catastrophe theory, 
instances of saddle connections.) 



B. Stokes transitions in the diffraction integrals 

In the neighbourhood of codimension two bifurcations 
one encounters tangent bifurcations in which satellites 
become ghosts, as was demonstrated in section [II. Sub- 



sequently the ghosts may undergo Stokes transitions with 
the central orbit. We analyze the transitions by investi- 
gating the expressions for the actions of the central orbit 
and the ghosts and using the condition (|5^). We will 
demonstrate that no connection can exist between the 
sheet of the central orbit and ghosts that have (almost) 
real J < and real actions. Moreover it is not difficult to 



realize that the condition (56) is on the other hand suf 



ficient for ghosts beyond tangent bifurcations at positive 
/. The reasoning is facilitated by the observation that 
the problem can be reduced in all cases to one dimension 
by considering the /-lines (or, to be precise, the complex 
/-planes) 4> mod tt = const that connect the central orbit 
and the ghost satellites radially. 

For m — 1 and m = 3 the situation is simple since there 
are no other orbits than the two satellites (which arc the 
ghosts in question) and the central orbit. The diffraction 
integral for to = 1 involves Pearcey's integral (and its 



derivatives) for which the Stokes transitions have been 
studied by Wright |^^. In the construction of the 
real solution has been placed in the origin, = po = Q. 
The satellites have coordinates p± — Q, 



I 'ia % o? le 



(57) 



(A tangent bifurcation is now encountered at £ = 
9a? /2>2h. Two orbits coalesce also at e = 0, but both 
remain real there due to the construction.) The orbit at 
the origin has the action S'o, while 



'5'± — So 



A l2 + 4'Z± 



(58) 



for the satellites. The Stokes transition takes place at 



(59) 



Figure displays the integration contour in the complex 
q'-plane for e — 3(3 -f \/3)/16, a = h = \ together with 
the equipotential lines of |exp[i$]| (or, equivalently, of 
Ini$). The plot demonstrates the well-known existence 
of the connection and is characteristic even for arbitrary 
sets of parameters that fulfill ( [59| ) since the shape of the 
contour is fully determined by the combination ehjc?: 
The contour expands linearly with a scaling of g, such 
that we can achieve, for instance, a — b, and does not 
change if ^^-^"^ is multiplied by a real constant, which 
allows to set b = 1. 

For even m the problem is mapped on the case to = 2. 
We already determined the locations ( pO| ) of the orbits 
for the Hamiltonian normal form H^^K By convention, 
coefficients changed sign in the definition of but the 
coordinates are not affected by this. The satellites have 
the actions 



la 2 a3 „, /la^ le 
^- = ^"+3^-27 6^±^n96^-35 



3/2 



(60) 



From the condition Rc S± = Sq one finds a Stokes transi- 
tion of complex ghosts for a = if eb > 0. (No transition 
is encountered at e = |^ or for a = 0, but eb < 0, since 
the radicand is positive then). As for to = 1,3 there 
is only one scale-invariant parameter combination eb/a^, 
but for each value there are now two variants of the con- 
tour depending on sign(a6). The Stokes transition at 
a — involves a ghost with a nonvanishing imaginary 
part of the action. Figure ^ shows how the integration 
contour changes in the complex q'-plane due to the tan- 
gent bifurcation and the Stokes transition. The plots 
proof the existence of a path that connects both sheets 
at constant Re$. 

Different situations with Re S± ~ So appear, however, 
1 ^ 

for £ = JIT- Stokes transition happens there because 
the satellite involved is real for ab < and a ghost with 
real action and real / < for ab> 0. Figure [lC| confirms 



13 



that indeed the contours are separated by another 'real 
ghost' (dashed contour). 

For TO = 5 the phase function on the p-axis is a polyno- 
mial of degree five. It appears on first sight that the three 
satellites could be arranged in such a way that the sheet 
of the ghost + is separated from the sheet of the central 
orbit by the remaining real satellite 1. One easily finds, 
however, that the sheets are separated for the particular 
phase function only if the real parts of the three roots p± , 
Pi of equation ( |2^ ) all have the same sign, which in turn 
can be ruled out by a careful inspection of equation (23): 



From Vieta's relations it follows that otherwise the coeffi- 
cient of the linear term ~ [p+P- + (p+ +P-)pi\ would not 
vanish. Note that the order I^/"^ in <i>*^^) indeed results in 
a nonvanishing coefficient but is considered, as usually, 
only as a perturbation and does not alter the situation 
qualitatively. The existence of a steepest-descent con- 
nection of the ghost and the central orbit at the Stokes 
transition is then guaranteed also for m — 5. The tran- 
sition takes place at the value oi t = 75e6^/(4a'^) that is 
the solution of 



640 



1520 i 925 



= 



(61) 



close to t — —97.6566. To derive this equation we first 
reduce the expression for Re S~Sq with help of the fixed- 
point equation 



(62) 



for the scaled variable r ~ p to the form Re [300r^ + t{8 — 
10r-|-75r^)] = 0. Introducing r = (x + iy) here and in the 
fixed-point equation and splitting the latter into real and 
imaginary part, one can solve then for y'^ = 8x/5 + 3.t^ 
and t = 15a;(2 + 5x)'^ and obtains the cubic equation 
8 -I- 60x + 125x^ + 75x^ — 0. Its roots give three values 
for t which all solve (p2|). Only one of the roots, however, 
fuffiUs = 8a;/5-f 3 > 0; this is the one that corresponds 
to the approximate value of t given above. 

For odd to > 7 there are four different satellites on each 
line. They are generically grouped in pairs that lie on 
opposite sides of the center, and the situation is similar to 
TO = 2 with slightly broken reflection symmetry. Indeed 
it follows from the derivation of the normal forms that the 
odd terms originating from higher-order perturbations 
are neglegible for small e and a; the results for m = 2 are 
then directly applicable to the present case. 



VII. CONCLUSIONS 

We studied bifurcations of codimension two in Hamil- 
tonian systems that are either autonomous and have two 
degrees of freedom or periodic with one degree of free- 
dom. The normal forms derived in section O and dis- 



cussed in section [II show that the typical sequences of 
codimension one bifurcations in the neighbourhood (in 
parameter space) of the bifurcation of codimension two 



consists of a period-TO bifurcation at a central orbit fol- 
lowed by a tangent bifurcation in which satellites become 
ghosts. 

Additional generic scenarios are encountered in the 
presence of symmetries []l6|-p0t. Isochronous pitchfork 
bifurcations are the most important addition of codimen- 
sion one in the case of time-reversal or reflection symme- 
tries; they will also show up in the neighbourhood of 
codimension two bifurcations in these systems. 

Only a small number of the bifurcations of codimension 
one and two correspond to (special cases of) a so-called 
elementary catastrophe due to Thom (see e.g. [ |37| , p8[ ). 
These appear in many different contexts and describe, for 
instance, bifurcations of codimension up to four in maps 
that are not restricted by area preservation. We use the 
usual names and symbols and further denote each Hamil- 
tonian bifurcation type by (toj,), where to is the multi- 
plicity and k = 1,2 the codimension. The fold A2 corre- 
sponds to the tangent bifurcation (li). The cusp A3 is 
(I2), and (2i) is a cusp with a reflection symmetry. (22) 
is a butterfly Ar, with reflection symmetry. The period- 
tripling bifurcation (3i) corresponds to a version of the 
elliptic umbilic DJ. All other normal forms describe 
catastrophes that would be of much higher codimension 
without area preservation. Especially for the cases to > 3 
one has to rely on higher-order perturbation theory. It 
implies that i) for a given codimension the class of bi- 
furcations in Hamiltonian systems is considerably larger 
and ii) although this can be circumvented by considering 
a normal form of much higher codimension from ordinary 
catastrophe theory, these normal forms have then again 
to be restricted: Points that correspond to the trajectory 
of one and the same orbit lie on the same height (energy 
or action). The classical perturbation theory takes care 
of this and in addition gives the right codimension. 

Collective contributions to semiclassical traces were de- 
rived that involve normal forms for a phase function $ 
and an amplitude function ^P. The expressions involve 
just as many coefficients as are determined by the ac- 
tions and stability properties of the bifurcating orbits, 
including a suppression of certain unwanted ghosts for 
TO > 4. The expressions constitute uniform approxima- 
tions: They are also valid far away from the bifurcation 
and asymptotically take the form of a sum of isolated 
contributions (^9|). The uniform approximations display 
Stokes transitions in which the ghost satellites interact 
once more with the central orbit and leave the steepest- 
descent integration contour. 

The validity of the approximations given here is lim- 
ited if additional orbits become important or 'unwanted' 
ghosts become real; bifurcations of even higher codimen- 
sion are then to be studied. The basic steps would be 
the same as in the present study: Derivation of Hamilto- 
nian normal forms that account for all bifurcating orbits; 
reduction of normal forms to get rid of non-bifurcating 
orbits and to account for independent stabilities and ac- 
tions. An important open question is concerned with 
the complexity of periodic-orbit clusters typically en- 
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countered in the quest of resolving spectra when one ap- 
proaches the semiclassical hmit: This requires knowledge 
of the dynamics up to the Heisenberg time ^ 1/h and 
involves a competition of increasing resolution in phase 
space and proliferation of periodic orbits. 

One might also be concerned about cascades, which 
are sequences of bifurcations of a certain orbit of pe- 
riod n at differing values of trM*^"' = 2cosci; (cf. section 
II B). The most prominent example is the basic build- 



ing block of period-doubling cascades: an orbit of period 
2n is born at an orbit of period n in a period-doubling 
bifurcation (trM*^^") = 2) and period-doubles itself at 
^^,j^{2n) _ _2_ ^ huge variety of cascades exists, how- 
ever, since bifurcations happen whenever the stability 
angle w is a rational multiple of 27r. Bifurcations in a 
cascade cannot be encountered simultaneously in param- 
eter space, since this would imply a singular change in 
the linearized map. For that reason cascades cannot be 
regarded as unfoldings of bifurcations of higher codimen- 
sion. The bifurcations in an unfolding show up simultane- 
ously in a given iteration of the map; the cascades involve 
bifurcations that appear in distinct iterations. One could 
study, for instance, those cascades that arise from the it- 
eration of the map generated by a normal form, and ask 
the question whether situations exist in which the orbits 
in the cascade must be treated collectively; it would be 
indeed nice to see that one can do without. An argument 
in favor of this expectation has been given in . 
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FIG. 1. Contour plots of the normal form i?'^' as the parameters are steered to cross two tangent bifurcations close to the 
bifurcation point of codimension two. Initially, only one orbit is present. Two new orbits are born in a tangent bifurcation. 

One of them approaches the first orbit, and both annihilate in an inverse tangent bifurcation. A similar scenario exists in which 
stable orbits are unstable and vice versa; it is obtained by reversing the sign of a. 




FIG. 2. The typical sequence for m = 2 of a tangent bifurcation of period-two satellites and a period-doubling bifurcation 

is illustrated by contour plots of H^'^\ As for m = 1 there exists a similar scenario for the opposite sign of g in which the 
stability of orbits is changed. The tangent bifurcation would not be encountered in real phase space if the satellites meet at a 
negative value of /. 




FIG. 3. The contour plots of the normal form H^^^ display a sequence of a tangent bifurcation of satellites and a pe- 
riod-tripling bifurcation. 




FIG. 4. The two sequences of contour plots of H^'^^ display a tangent bifurcation of satellites followed by a period-quadrupling 
bifurcation. The latter is encountered in its island-chain version above; below we have the touch-and-go scenario. 
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FIG. 6. Varying a parameter close to the codimension two point with m = 6, one might observe two pairs of satelHtes being 
born in a tangent bifurcation at positive I, as displayed here; the satellites closer to the center disappear in a subsequent 
period-6 bifurcation. The island chain that is left over here could be also steered to the center by letting a, and then once 
more s change its sign. 




FIG. 7. For m > 7 the resonant </)-dependence of the normal form is weaJi, and tangent bifurcations of satellite pairs happen 
at almost identical values. Subsequently the inner island chain collapses onto the center and disappears. The remaining chain 
might follow, as explained for m = 6. 




FIG. 8. Path of steepest descent Re 'I>''^'' (g, 0) = So (thick line). The parameters are chosen to fulfill the condition for 
a Stokes transition. The transition indeed takes place since the contour connects the subdominant ghost with the dominant 
central orbit. The thin lines are the equipotential lines of | exp[i$]| (or Im$). 



17 



FIG. 9. Sequence of steepest-descent contours for {q, 0) displaying a tangent bifurcation and a Stokes transition. Again 
there is a connection of the ghost satellite and the central orbit as the condition for a Stokes transition is fulfilled. Dashed 
lines indicate steepest-descent paths that are not needed to connect the integration boundaries. A trick can be played with 
these pictures to envision the situation for the tangent bifurcation at negative I: The plots are rotated by 90 degrees (which 
corresponds to inverting the sign of a) , and the contour is picked that originally connects ±ioo (see also Fig. uXm . 




FIG. 10. No Stokes transition happens for ghosts with real 7 < in the representative case m = 2. The steepest-descent 
contours from the ghost and the central orbit (full lines) are separated by the (dashed) contour of another ghost at real, but 
negative /. Only the contour that visits the central orbit is needed to connect the integration boundaries. 
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